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ON  EVALUATION  OF  SOLUTIONS  TO  EQUATIONS  OF  WAVE 
PROPAGATION  IN  A  LAYERED  HALFSPACE 


ABSTRACT 

To  remove  certain  discrepancies  in  solutions  of  problems  of  wave 
propagation  in  layered  media,  especially  in  the  case  of  a  two-layered 
liquid  halfspace  first  discussed  by  Pekeris,  the  precise  meaning  of  repre¬ 
sentations  used  is  reconsidered.  The  existence  of  real  poles  on  the  path 
of  integration  used  by  Lamb  and  others  for  real  CJ  leads  to  the  addition  of 
corresponding  residues  to  the  principal  value  of  the  integral.  They  repre¬ 
sent  the  free  waves  usually  superposed  on  the  solutions  to  yield  divergent 
waves.  The  question  of  real  and  complex  poles  is  clarified  for  the  Pekeris 
problem  and  the  new  form  of  solution  is  given  in  terms  of  a  branch  line 
integral  and  residues  at  all  these  poles.  The  Lamb  condition  of  vanishing 


potential  at  infinity  is  fulfilled. 


2. 


1.  The  solutions  of  problems  on  wave  propagation  in  layered  media 
by  the  well  known  methods  of  Lamb  (1904)  and  Sommerfeld  (1909,  1949)  are 
usually  obtained  first  in  the  form  of  certain  integrals.  These  integrals  can 
be  improper  and,  therefore,  use  is  to  be  made  of  their  principal  values  in 
the  Cauchy  sense,  as  has  been  done  by  Lamb  in  the  case  of  a  homogeneous 
half  space.  It  was  pointed  out  by  Kupradze  and  Sobolev  (1930)  that  the  pro¬ 
cedure  used  by  Lamb  is  equivalent  to  the  choice  of  a  definite  value  of  a  di¬ 
vergent  integral.  Further  transformations  of  a  solution  of  a  wave  problem 
lead  to  other  difficulties  of  mathematical  character  which  are  not  avoided 
in  certain  modifications  of  Lamb's  method.  For  example,  Pekeris  (1948) 
derived  expressions  for  potentials  representing  wave  propagation  in  two¬ 
layered  liquid  halfspace  in  terms  of  an  infinite  number  of  normal  modes  and 
of  a  branch  line  integral.  Recently  Stone  (1953)  pointed  out  that  Pekeris 
did  not  include  in  his  solution  an  infinite  number  of  complex  poles.-  Since 
these  poles  are  located  in  the  fourth  quadrant  they  yield  in  his  solution 
terms  exponentially  increasing  with  depth.  It  is  obvious  that  the  existence 
of  these  terms  presents  a  great  difficulty  to  the  theory  and  Stone  derived 
another  form  of  solution  of  the  two-layer  problem. 

The  existence  of  roots  of  the  period  equation,  i.  e.  poles  of  integ¬ 
rands,  was  investigated  by  Sherman  (1945)  in  more  detail  for  the  case  of  a 
liquid  layer  over  a  semi-infinite  solid.  He  proved  that  there  are  an  infinite 
number  of  simple  complex  roots  of  the  period  equation,  besides  a  finite 
number  of  real  roots,  and  pointed  out  that  the  point  at  infinity  is  an  essen- 


3. 


tial  singularity  of  the  integrands  obtained  in  this  problem. 

The  fact  that  there  are  certain  discrepancies  in  representations  of 
solutions,  one  of  them  being  just  mentioned,  requires  special  attention. 

We  have,  therefore,  to  make  a  revision  of  the  approach  to  certain  points 
of  the  theory.  The  solution  of  the  problem  of  wave  propagation  in  a  two- 
layered  liquid  halfspace  will  be  used  as  an  illustration. 


2.  We  begin  with  some  remarks  concerning  the  expressions  repre¬ 
senting  potentials  due  to  a  harmonic  point  source  in  a  liquid  layer  overlying  a 
semi-infinite  liquid.  We  write  the  expressions  obtained  by  Pekeris  (1948) 
in  the  form 


a inh  y'z  coi/t  v '(H-h )  p-, si fii  v'(H-h) 


V' 


A 


azs/j  y  '(z-  H)  -  y  '(z  -H) 

A 


z<k  (i) 

z>k  (2) 

(3) 


where 


V'H  +  6 * 


SinJi  y' H 


Pi 

H  is  the  thickness  of  the  layer  and  &  =  is  the  ratio  of  densities.  As 

ft. 

usual  r  and  z  are  cylindrical  coordinates  and  k  the  wave  number.  In  the 
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potential  ^  for  the  semi-infinite  liquid  substratum  ( //=  Z  ~  )  we 

have  a  factor  exp  l-y"  (z-H)]  .  In  order  to  obtain  functions  which 

will  not  diverge  as  z  -►C®  it  is  assumed  that  for  a  real  k  >  —  and 

k  >—  respectively 


*  /<■*-  s 


2. 

£ 


(5) 


are  real  and  positive.  In  the  problem  we  are  considering  the  source  is 
located  at  the  point  (0,  0,  h)  of  the  layer. 

The  expressions  (1-3)  have  characteristics  typical  of  solutions  of 
similar  problems.  The  Bessel  function  J  (kr)  in  (1)  is  usually  replaced  by 


,0) 


Hankel  functions:  Z  Jc  («r)  =  H"(«r)  +  H(*W)  and  in  further  trans¬ 


formations  either  H 


(i) 


or  H 


o 

(2) 


cancel  out.  The  reason  for  the  choice  of 


o  o 

one  of  these  functions  is  immediately  seen  from  their  asymptotic  form 


Z  12  < 
nurj  ® 


provided 


ft)'  » f  Z  -tVr 


/rVr/ 


provided  —  Z.H  +  ZS ~  ^  tr -2 s 


S  being  any  positive  acute  angle.  If  one  makes  use  of  the  time  factor  in 


the  form  e***^  the  second  Hankel  function  must  be  taken  in  the  final  expres¬ 


sion. 


To  these  preliminary  remarks  we  add  the  following  ones.  Since 


the  period  equation 


A  =  c*/,  y'H+  gy  ^  o 


(6) 
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yields  in  general  a  set  of  roots  the  expressions  (1-3)  can  be  improper  in¬ 
tegrals  because  of  the  existence  of  real  poles  of  the  integrands.  Similar 
expressions  are  also  improper  integrals  because  of  infinite  limits.  There¬ 
fore,  we  will  use  the  notation  for  their  principal  values  in  further 

transformations. 

In  (1-3)  the  summation  is  over  the  range  of  positive  values  of  the 
real  variable  k  in  which,  besides  the  poles,  there  are  branch  points  of 
two  radicals  (5).  We  see,  however,  that  the  change  in  sign  of  Yf  does 
not  affect  the  expressions  (1-3)  and,  therefore,  the  corresponding  branch 
line  integral  vanishes  as  was  found  by  Pekeris  (1.  c).  A  generalization  of 
this  important  property  of  solutions  for  an  n-layered  halfspace  was  given 
by  Jardetzky  (1953). 

According  to  Sherman's  analysis  mentioned  above  we  can  expect 
that  the  period  equation  (6)  has  a  finite  number  of  real  roots  in  the  interval 

/CO  CaJ_ 

(  <*2  ?  c*., 

since  this  is  a  more  interesting  case  (a  "£a&  bottom"),  and  well  known 
simple  considerations  determine  the  number  and  location  of  real  roots  in 
it.  As  to  the  complex  roots,  Sherman's  proof  is  based  on  the  rule  concerning 
the  number  of  zeros  and  poles  of  a  function  of  complex  variable.  We  will 
not  yet  specify  their  number,  mentioning  first  the  conditions  under  which 
the  evaluation  of  integrals  of  the  type 


)  .  We  usually  consider  this  interval  assuming  that 


(?) 
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by  contour  integration  are  proved  to  he.correct.  The  integrals  in  (1-3) 
have  this  form  and  the  function  Q  is  also  a  function  of  the  set  of  variables: 
r,  z,  ,  A, 

It  is  easy  to  see  that  the  integrands  in  (1-3)  can  also  be  considered 
2 

as  functions  of  k  ,  i.  e.  can  be  replaced  by  an  integral  of  the  form 


= /  Q(*)M 


(8) 


-  CO 


Now  it  has  been  proved  that  the  principal  value  of  such  an  integral 
can  be  computed  by  Cauchy  theorem  and  is  given  by  the  formula  (for  k  real) 


=  2ri(ZTt+jZT]0) 

~oo 


(9) 


where 


means  the  sum  of  residues  of  the  poles  on  the  real  axis. 


2L  R  is  the  sum  of  other  residues,  if  the  function  Q(k)  satisfies  certain 
conditions  (Whittaker  and  Watson,  1947,  p.  111).  These  conditions  are  as 
follows.  Let  us  consider  a  complex  variable  £  *  tC  ■+  ^  and  a  function 
Q(  Zy  )  such  that:  l)  it  is  analytic  when  the  imaginary  part  of  Z»  is  positive 
or  zero  except  at  a  finite  number  of  poles;  2)  it  has  a  finite  number  of  simple 
poles  on  the  real  axis;  3)  as  /£/-><*>,  £  Q(&)-+  0  uniformly  for  all 

values  of  arg  £  such  that  0 ssr  =  7T  and  4)  when  k  is  real 

kQ(k)  — >■  0  as  k-^  ±  CO  in  such  a  way  that  J  and  f  both 

converge.  If  there  are  trigonometric  factors  in  Q  the  integrands  can  be 
subjected  to  less  stringent  conditions.  The  condition  that  Q(  £  )  is  a  single- 
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valued  function  must  be  attended  too  (Whittaker  and  Watson  1.  c.  p.  113). 

These  conditions  are  given  in  a  form  appropriate  to  the  application  of  the 
Cauchy  theorem  in  the  upper  halfplane.  They  can  be  easily  modified  when 
the  integration  occurs  over  the  lower  halfplane. 

In  this  connection  a  remark  has  to  be  made  concerning  the  use  of 
the  complex  £  plane.  Because  of  the  existence  of  branch  points  certain 
cuts  must  be  used  to  determine  the  corresponding  Riemann  surface.  Then, 
for  example,  the  admissible  first  and  fourth  quadrants  can  belong  to  differ¬ 
ent  sheets  or  to  the  same  primary  sheet  for  a  multi-valued  function.  This 
fact  can  easily  produce  further  confusion,  if  the  integrand  has  singularities 
and  if  these  points  belong  to  different  sheets  of  a  Riemann  surface.  By  the 
choice  of  a  cut  one  can  omit  or  introduce  some  singularities  which  are  not 
appropriate  to  the  solution. 

It  has  been  proved  (Jardetzky,  1953)  that  in  a  problem  of  wave 
propagation  in  an  n-layered  medium  all  the  integrands  (like  1-3)  are  even 
functions  of  all  (and  for  shear  waves)  except  for  those  exponents 

which  refer  to  an  infinite  part  of  the  medium  (the  last  two,  Vn  and  Vn^ 
in  case  of  a  halfspace).  Therefore  all  cuts  but  four  referred  to  the  branch 
points  ±  ~  and  ±  can  be  omitted.  A  s  to  these  four  cuts  they 

<*n  Pn 

should  be  drawn  always  in  such  a  way  that  no  roots  of  the  period  equation 
which  violate  Lamb's  condition  are  on  the  chosen  part  of  the  Riemann  surface. 
This  condition  in  a  general  case  is  that  no  negative  real  parts  of  the  coef¬ 
ficients  or  of  a  positive  z  in  exponential  factors  exp(-Vz)  of  in- 
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tegrands  are  permitted.  All  these  four  branch  points  appear  of  course  in 
the  case  of  a  homogeneous  solid  halfspace  (Lapwood,  1949)  and  were  con¬ 
sidered  also  in  the  problem  of  a  liquid  layer  over  the  solid  halfspace 
(Sherman,  1945,  compare  Newlands,  1952). 

3.  Now  writing  the  integral  (l)  or  a  similar  one  in  the  form 

<f  =  («r) F(v'<  * ') < i o) 

0 

we  have  to  see  the  precise  meaning  given  to  it.  There  is  an  essential 
difference  in  problems  considered  by  Lamb  (1904)  and  Sommerfeld  (1909). 
The  latter,  investigating  the  propagation  of  electromagnetic  waves,  ob¬ 
tained  complex  branch  points  and  a  complex  pole  so  that  in  his  solution 
there  are  no  singular  points  of  the  integrand  in  (10)  on  the  real  axis. 

Thus  no  special  definition  of  the  integrals  used  was  needed.  A  similar 
situation  is  found  in  application  of  Sommerfeld’ s  method  by  Lapwood 
(1949)  to  elastic  waves  when  a  complex  U)  and  operational  analysis 
are  used  instead  of  a  real  frequency  and  Fourier  analysis.  For  real  CJ  , 
however,  beginning  with  Lamb's  problem  we  find  in  problems  concern¬ 
ing  wave  propagation  in  layered  media  singularities  on  the  real  k  axis. 

A  precise  definition  of  expression  (10)  is,  therefore,  required  and  this 
can  be  obtained  from  the  discussion  given  by  Whittaker  and  Watson  (1.  c.  K 
The  integration  in  (10)  has  to  be  performed  along  the  positive 
part  of  the  real  axis  between  the  poles  and  along  infinitely  small  half- 
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circles  with  centers  at  these  points.  Thus  if  X+n  are  the  poles 


Y,ur\f+r\f+...\ 

r  Vo  r,  Vi  } 


provided  the  radii  V  tend  to  zero.  Now  the  sum  of  odd  terms  is  the 

tv 

principal  value  of  the  integral.  The  terms  taken  along  the  half-circles 

^  yield  the  residues  corresponding  to  real  poles  which  are  obtained 

T  p  1  ^ 

by  the  transformation  &  -  **  ^ 


and  lim  lQ(Z)(Z-x,j]  =  Res.  Thus 

2J„(*r) F(y[  -  riZ *R,  ( 1 1) 

0 

where  the  sum  is  taken  over  a  finite  number  of  real  poles.  Equation  (11) 
yields  the  explanation  of  the  fac£  that  in  order  to  obtain  a  train  of  waves 
diverging  from  the  source  Lamb  had  to  add  "free  Rayleigh  waves"  to  that 
part  of  his  solution  which  is  represented  by  principal  parts  of  integrals. 
This  train  of  "free"  waves  is  exactly  given  by  the  last  terms  in  (11). 

An  attempt  has  been  made  by  Lapwood  to  explain  the  addition  of 
free  waves  as  due  to  the  special  choice  of  the  Riemann  surface  made  by 
Lamb.  It  is  obvious  that  choosing  the  path  of  integration  as  a  real  axis 
under  other  conditions  we  will  obtain  different  definitions  of  the  integral 
(10).  Especially,  using  different  rules  for  construction  of  a  Riemann 
surface,  we  will  have  to  vary  the  contours.  However,  the  relationship 
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(11)  is  due  to  the  existence  of  real  roots  of  the  period  equation  and  holds 
for  a  single- valued  function  as  well  as  for  a  multivalued  integrand  for 
which  the  use  of  a  Riemann  surface  is  necessary  for  contour  integration. 
The  superposition  of  free  waves  appears  now  as  a  natural  part  of  the 
precise  definition  of  an  integral  representing  the  solution  of  a  problem  of 
wave  propagation  and  loses  the  character  of  an  arbitrary  transformation. 

4.  The  second  point  in  the  definition  of  (10)  which  requires  more 
precision  is  to  obtain  a  uniform  function  as  an  integrand,  as  mentioned 
before.  In  the  case  of  the  integral  (10)  we  have  a  two -leaved  Riemann 
surface  because  of  two  values  of  the  radical  V11  .  It  is  important  to  indi¬ 
cate  all  conditions  under  which  the  choice  of  such  a  surface  is  made,  es¬ 
pecially  those  concerning  an  appropriate  cut.  Sherman  and  others,  for 

example,  used  as  a  cut  the  part  of  the  real  axis  between  the  branch  points 
,  Co  \ 

±  =±  ■«£-  (see  equation  3).  Pekeris  made  use  of  rays  starting  at 

the  branch  points  and  k^  and  parallel  to  the  negative  imagin¬ 

ary  axis.  Such  cuts  seem  to  produce  no  confusion  in  Lamb’s  problem  but 
they  were  also  used  in  many  other  investigations.  This  choice,  accord¬ 
ing  to  Stone  (1953),  is  responsible  for  inclusion  of  terms  having  no 
physical  meaning.  In  discussions  of  wave  propagation  in  different  cases 
either  upper  or  lower  and  also  the  left-hand  or  right-hand  halfplanes  were 
used  as  a  part  of  the  admissible  sheet  of  a  Riemann  surface.  We  will 
take  the  latter.  To  get  more  insight  into  all  these  questions  it  is  easier 
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to  begin  the  discussion  for  complex  <o  as  it  was  done  by  Lapwood. 

Since  in  potentials  for  the  infinite  part  of  a  medium  we  have  a 
factor  exp  (  -  Y'1  z'j  Rev"  cannot  be  negative  for  positive  z  as  has  been 


assumed  above.  Now  putting 

y's  a-h  ?  ~  X  +  , 


CO 


*°I~LP 


we  can  find  the  curve  determined  by  the  condition 

7tc  y"~  a-  -  0 


By  (5) 


-  /*=  kz-  rV  J? Vfi  - 


1*'-  P*'-  2<]Pl 


(12) 


(13) 


(14) 


Thus  the  line  in  question  is  a  part  of  the  hyperbola  (Fig.-  la). 

(is) 


CJ 


It  starts  at  the  branch  point  £  ^  in  the  fourth  quadrant,  if  we 

2 

restrict  ourselves  now  by  the  case  q>0,  p  ^  0.  This  follows  from  the 
second  condition  given  by  (12),  namely, 


yir*<  S-=JSi 


(16) 


We  can  easily  see  that  the  second  part  of  the  hyperbola  (15),  that  on  the 


other  side  of  the  branch  point,  corresponds  to  the  conditions 

J/n  r'=  £  =  0 


(17) 


and 


A"*-  Zi  > 


<x. 


(18) 
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At  the  branch  point  k  the  real  and  imaginary  parts  of  V^both 
vanish  and  a  similar  distribution  around  the  branch  point  -k^  occurs 
on  the  second  branch  of  the  hyperbola  (15)  in  the  second  quadrant.  A 
change  in  sign  of  p  or  q  will  displace  these  curves  into  first  and  third 
quadrant  respectively. 

There  are  two  extreme  cases:  CO  =q  is  real  and  U>  =-ip  is 
purely  imaginary,  in  which  (15)  becomes 


xv=  0 


(19) 


For  CJ  real  we  obtain  by  ( 1 5  -  1 6) :  1 )  T?-0  ,  k  *  <  -^3 


a, 


*  V.  ** 


i.  e.  a  part  of  the  real  axis  between  the  branch  points,  or  2)  k=0,  -“T  < 

i.  e.  the  imaginary  axis.  Taken  separately  this  axis  does  not  represent 
a  cut  since  it  does  not  pass  through  a  branch  point.  It  can  be,  however, 
combined  with  the  line  BOA  (Fig*  lb)  to  form  two  cuts  AOD  and  BOC 
which  obviously  represent  the  degenerated  branches  of  the  hyperbola  (15) 
if  p->  0  (Lapwood,  1949). 

If  we  take  CO  =  -ip,  p>0,  i.  e.  the  time  factor  e1****  =  e^  used 
first  by  Cagniard  ( 1939)  we  obtain:  l)  k=0,  T*  >  p*"  i.  e.  cuts  along  the 
imaginary  axis  or  2)  T  =0,  k^  ^  -p^,  i.  e.  conditions  which  cannot  be 


satisfied. 


We  will  now  consider  the  case  of  Cj  real  in  more  detail.  The 
line s  AO  D  and  BOC  are  determined  by  the  condition  (13)  and  Re  V  can 

either  change  sign  or  remain  positive  or  negative  on  both  sides.  Now  for 

// 

a  positive  z  we  have  to  associate  a  positive  Re  ^  with  poiifts  to  the  right 
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of  the  point  A  on  the  real  axis.  This  first  condition  will  determine  the 
admissible  sheet  of  the  Riemann  surface  for  this  case.  Now  it  can  be 
shown  that  under  appropriate  conditions  the  right-hand  halfplane  of  £ 
belongs  to  it.  Let  us  put 


it 


Then  by  (3) 


(20) 


(21) 


and  in  the  neighborhood  of  A  on  neglecting  T  we  obtain 

Hence 

V'=  -  e  <r|  +  r; 


If  we  assume  that  in  the  right-hand  halfplane  the  principal  value 
of  the  angle  ^  is  used,  i.  e.  -  7T<  $  £  7T  it  is  the  first  definition  (22) 
of  the  radical  which  yields  a  positive  Re  V*  in  this  halfplane.  Obviously 
Im  V  is  then  positive  in  the  first  quadrant  and  negative  in  the  fourth. 

To  insure  the  vanishing  of  exp(—  V"  z)  at  z  ->nTO  for  z<  0  we  should  use 
another  convention  of  signs  (e.  g.  the  second  definition  (22))  or  change 
the  sign  in  the  exponent.  In  case  of  a  complex  co  Im V,f  is  positive  on 
one  side  of  the  hyperbola  and  negative  on  the  other  (Fig.  la).  On  a  part 
of  the  hyperbola  or  on  the  part  Ak  of  the  real  axis  (Fig.  lb)  where  Im 
Y1'  is  continuous.  It  is  discontinuous  along  OA. 
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If  BOC  or  BOA  is  the  cut  we  put  for  the  neighborhood  of  point  B 

and  vary  the  angle  ^  in  the  interval  o  ±  2tr  in  order  not  to  cross 


the  cut.  Then 


•  r - i±pi  ,,  r - „ 


rf 


(f 


(22)) 


Re  Y,  is  now  negative  if  the  first  definition  is  taken. 

For  the  further  transformation  of  the  integral  (10)  we  need  only 

the  right-hand  halfplane  and  it  is  important  to  note  that  under  the  condi- 

_  # 
tion  that  the  angle  v  has  its  principal  value  the  definition  of  Y  along 

OA  is  not  changed.  The  change  in  sign  of  the  radical  occurs  only  for  its 

imaginary  part,  the  Re  v"  being  zero  at  every  point  on  OA  and  positive 

on  both  sides  of  it.  At  A  tjeV  =0  and  ImV  =o. 

Now  using  this  definition  of  Y  the  integrand  in  (10)  becomes  in 

the  right-hand  halfplane  a  single -valued  function  of  Z.  as  is  required  for 

the  contour  integration.  If  the  integrand  is  a  continuous  function  of 

ft  .  h 

it  is  also  a  continuous  function  of  Re  Y  but  a  discontinuous  one  of  Im  Y 
along  the  cut  since  on  different  sides  of  it  the  radical  which  becomes 
equal  to  its  imaginary  part  has  opposite  signs.  When  real  <0  is  used 
the  part  OA  of  the  cut  coincides  with  a  part  of  the  integration  path.  To 
give  a  precise  meaning  to  the  integral  we  have  to  take,  for  example,  the 
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value  of  V  on  the  upper  side  of  the  cut.  Given  by  Lamb  as  a  part  of  a 
convention  on  signs  of  radicals  this  seems  to  be  arbitrary  because  such 
a  formal  choice  of  signs  does  not  show  the  true  reason  for  it.  The  exist¬ 
ence  of  the  discontinuity  along  the  cut  explains  now  the  appearance  of  the 
branch  line  integral. 

If  we  put  V  =a  the  real  and  imaginary  parts  can  be  computed 

for  every  given  point  Z»  from  (22).  Then,  since  a^  and  -b^  are  roots 
of  a  certain  quadratic  equation,  we  obtain,  for  example,  a=±Y^  .  The 

convention  to  be  attended  is  that  Re  V  >0  in  the  first  and  fourth  quadrant, 
i.  e.  a=  JiT, 


5.  We  mentioned  the  fact  that  in  case  of  a  liquid  over  a  solid 
half  space  there  are  an  infinite  number  of  complex  roots  of  the  period 
equation.  Since  this  case  includes  the  problem  of  a  liquid  half  space  we 
can  assume  that  Sherman' s  proof  also  shows  the  existence  of  such  roots 
in  our  problem.  To  locate  them  we  can  use  a  simple  transformation  as 
follows : 

Writing  exponential  functions  in  (6)  we  obtain  the  period  equation 
in  the  form 


or 

zy'H  g-y"-  y 1 

€ 


cv  %  Y  ' 


(24) 
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This  is  equivalent  to 


Py'H-  4  -  n-2rft*  0 

n  ™  gy»+yl 


(25) 


Every  root  of  (25)  satisfies  (23)  and  vice  versa.  Substituting  (5)  and  a  com¬ 


plex  Z»  for  k  we  obtain 

6"' 


r - —  a  6'/ 

zh  /r-<(  -  ^ 


- 


<? 


-+'yJzz-Kz 

» v  s  V/ 


’ —  /?.  S Ml  =:  (? 


(25') 


We  will  assume  now  that  lzl^>  1  and  put  £  ,  writing  the 


lefthand  member  of  (25')  in  the  form 

gy- a  - 


*<<*  "\2 


n.  2ii 


(26) 


Expanding  this  expression  for  large  values  of  Izl  we  can  restrict  our¬ 


selves  by  considering  the  terms 


^  r[t  ~  Z*  ”V2?  -  n^7r‘ 


(27) 


An  even  more  rough  approximation  leads  us  to  the  conclusion  that  the  roots 


are  in  the  fourth  and  second  quadrant  respectively.  If  we  put,  namely, 


G*0«-?rfZ, 


-t  £  -  n.  Sri  «  0 

’*■  6 


(28) 


we  see  that  there  is  a  root 


TTfl  . 

TL 


(29) 
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provided  n  is  large  enough  to  make  the  modulus  of  a  large  number. 

This  was  the  condition  under  which  the  approximation  has  been  made. 

Since  the  definition  (5)  has  been  used  and  it  holds  for  the  fourth 
quadrant  of  the  admissible  sheet  of  the  Riemann  surface,  there  are  complex 
roots  on  it.  It  should  be  noted  that  there  are  no  roots  of  the  period  equation 
(6)  on  the  imaginary  axis.  If  we  put  this  equation  takes  the  form 


ck?S  at  H  lS  ^  ~  0 

where  d(  =  yjm JC*{  ^  yjtn7--*-  .  It  has  no  real  roots  m. 

The  change  in  sign  of  ,  i.  e.  the  use  of  the  second  definition 

of  the  radical  (22)  will  affect  the  form  of  the  period  equation  (6).  Since 

the  sign  of  the  second  term  of  this  equation  will  be  reversed  the  term  In  in 

(25)  will  be  positive.  Then  it  may  be  easily  seen  that  the  roots  of  the  new 

period  equation  are  located  in  the  first  and  third  quadrant  respectively. 

6.  The  discussion  of  the  integrand  in  (10)  has  shown  that  there 

are  three  basic  characteristics  in  case  of  CO  real:  l)  the  finite  number 

of  real  toots  in  the  interval  (k  ,  k<x/  );  Z)  infinite  number  of  complex 

roots  in  the  fourth  quadrant;  3)  the  part  of  the  integration  path  from  0  to  k 

0\Z 

is  a  part  of  the  branch  line.  On  its  two  sides  V  ^  has  purely  imaginary 
values  with  opposite  signs.  To  give  a  precise  meaning  of  the  integral  (1X1) 
we  must  use  one  of  them,  say  the  positive  one  in  (10)  and  (11). 

If  we  now  substitute  2 JQ(kr)=H0^|kr)f  HQ(^)(kr)  in  the  principal  value 
we  obtain 


(30) 
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In  the  first  integral  we  can  use  the  contour  enclosing  the  first 

quadrant  where  there  are  no  other  roots  than  those  on  the  real  axis 

(Fig.  2).  The  Cauchy  theorem  (see  equation  9)  yields 

i  too 


I = ^(WFCv',  K )  i  di 


(31) 


where  are  residues  at  real  poles  when  one  of  the  factors  is  the  first 
Hankel  function,  £  -Vi  f  and  the  subscript  +  shows  that  along  the  positive 
imaginary  axis  we  have  y  ff=s  Jm  y  "  >  0. 

To  transform  the  integral  I  in  a  similar  way  we  must  take  into 

L 

// 

account  the  fact  that  below  the  cut  ImV  is  negative.  Therefore  we  will 
add  and  subtract  J  X](«r)F(Y\  vl)KcU  , 


write 


1 ;  -2  V  )  - F(v ;  vlpdK 


(32) 


and  apply  the  countour  integration  to  the  sum  I1.  Then 
r  -too 

r-f  H*\Zr) F(y [vl)Zdl~ 2niZ T^-viZ 


(33) 


where  jL'R  is  the  sum  of  residues  for  real  roots  and  the  function  H  ^  is 

o 

one  of  the  factors.  R  is  a  residue  for  one  of  the  complex  poles  (  2*^  ) 

in  the  fourth  quadrant.  Taking  into  account  the  relationship 


Hi’ (Z)  - 

the  last  term  in  (31)  is  (  £ " ti  ) 

[  Tffanl <  )^r  = f  ’  Tfrctri)  )  tdr 


(34) 
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Now  adding  (31),  (34)  we  obtain 


Since  in  (ll)  the  residues  are  computed  for  the  factor  ZJQ  we  have 


2Z? 


(36) 


the  first  sum  in  the  right  hand  member  being  taken  for  the  Hankel  function 

(i)  fVi  q  o  q  r*  rvr\  rl  T-T  (2) 


H  '  ,  the  second  for 

o  o 


Then  by  (8),  (9),  (35)  and  (36) 
oo 


cn  =  Zf  Jo(^r)f(y'n  v")kc(k=  -2/T'Z 3*iZ ^  + 

*  0 

J ffir)[F(y>y+)  -F vf  rl'jjz  dl  +f  ffW(*r)jF(vtf)-F(vy^K<lK. 


-001 


(37) 


The  residues  Rq  and  R  have  to  be  computed  for  the  integrand 

which  for  the  real  poles  takes  the  form  /f’MRy'SJ* 
and  r"=^K"  .  The  sum  of  the  last  two  terms  is  the  branch  line  integ¬ 
ral  taken  along  the  cut  -mioo}  OA  .  a  represents  superposition  of  a 
continuous  wave  spectrum.  It  is  well  known  that  when  such  integrals  are 


evaluated  by  different  methods  they  are  approximated  by  a  certain  number 
of  waves  in  a  different  way.  As  to  the  function  F  its  form  is  seen  in 
three  expressions  (1-3)  which  hold  for  wave  propagation  in  a  two-layered 
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liquid  half  space.  The  integrands  thus  have  the  form  /£  ^  (%)  /  1 

According  to  the  well  known  formula  the  residue  at  a  simple  pole  (k  )  is 


In  (1-3)  -A  and  <  is  its  derrcitive  with  respect  to  Z  , 

being  the  real  roots  of  the  period  equation  ^(k^^O. 

The  solution  of  this  problem,  represented  by  three  equations  of 

the  type  (37),  differs  from  those  given  in  investigations  of  Pekeris  and 

Stone  quoted  before.  As  mentioned  above,  Stone  pointed  out  that  had 

Pekeris  introduced  the  second  sum  of  residues  in  (37)  which  correspond 

/( 

to  complex  roots,  they  should  be  associated  with  a  negative  Re  Y  because 
of  his  choice  of  the  cut,  and  therefore  would  violate  the  condition  at 
infinity  (  2.  00  ).  Bestdes  these  terms  there  are  in  (37)  only  a  finite 

number  of  residues  RQ  corresponding  to  real  roots.  The  branch  line 
tmrgTgj  in  (37)  differs  also  from  the  straight  line  from  k  to  k^  -  lea 

used  by  Pekeris,  and  also  from  that  given  by  Stone.  The  latter  obtained 


in  his  final  expression  the  integral  from  0  to  k  only  since  he  assumed 

that  Pekeris*  proof  of  the  fact  that  terms  like  the  last  in  (31)  and  the  first 

in  (33)  cancel  out  holds  also  for  his  new  choice  of  the  cuts.  For  the  cuts 

shown  in  Figure  lb,  however,  the  values  of  V JVn  V ^  Qn  the  imaginary 

axis  differ  in  sign  on  its  positive  and  negative  parts  as  we  have  seen  before 

and  indicated  by  the  subscripts  +  and  .  Comparing  the  distribution  of  Re  y  " 
n 

and  Im  V  given  in  Figures  la  and  2,  it  may  be  easily  seen  that  the  Pekeris 
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conclusion  about  the  vanishing  of  integrals  is  correct  since  it  corresponds 
to  the  case  la.  It  does  not  hold  for  case  2  considered  by  Stone.  There¬ 
fore  in  our  case  the  part  of  the  branch  line  integral  from  —  too  to  0  does 
not  vanish.  In  Stone's  solution  where  a  finite  number  of  normal  modes  is 
obtained,  the  sum  of  residues  such  as  the  second  in  (37)  is  omitted,  but 
using  Sommerfeld's  cuts  we  have  a  positive  ReY^  in  the  whole  right- 
hand  halfplane  which  is  associated  with  the  first  definition  (22)  of  the 
radical.  We  have  seen  that  this  halfplane  belongs  to  the  admissible  sheet 
of  the  Riemann  surface.  The  roots  on  it  do  not  violate  the  condition  at 
infinity. 

The  new  form  of  solution  of  the  problem,  i.  e.  the  equations  (37) 
representing  potentials  of  wave  propagation  in  a  two-layered  liquid  half¬ 
space,  obviously  does  not  affect  the  conclusions  for  the  field  involving 
phase  and  group  velocities  reached  by  Pekeris  and  others.  These  applica¬ 
tions  of  the  theory  are  based  on  consideration  of  the  period  equation. 

We  discussed  here  only  the  case  of  one  radical  which  makes  the 
integrands  double  valued  functions.  A  generalization  for  multi-valued 
functions  requires  no  essential  changes. 

It  is  a  pleasant  duty  of  the  author  to  express  his  hearty  gratitude 
to  Dr.  Frank  Press  for  very  useful  discussions  of  the  subject  of  this  paper. 


References 


Cagniard,  L.  Reflexion  et  refraction  des  ondes  seismiques  progressives 

(ThEse);  Paris,  Gauthier  Villars,  1939 

Jardetzky,  W.  S.  Period  equation  for  an  n-layered  halfspace  and  some 

related  questions;  Tech.  Rep.  on  Seismology,  Columbia 
Univ.  No.  29,  1953 

Kupradze,  V.  and  Sobolev,  S.  Sur  la  propagation  des  ondes  elastiques  a 

la  surface  de  separation  de  deux  milieux  ayant  des  propriEte's 
Elastiques  diffe'rentes ;  Pub.  Inst.  Seism.  Acad.  Sci.  USSR 
No.  10,  1930 

Lamb,  H.  On  the  propagation  of  tremors  over  the  surface  of  an  elastic 

solid;  Phil.  Trans.  R.S.A.  203,  1-42,  1904 

Lapwood,  E.  R.  The  disturbance  due  to  a  line  source  in  a  semi-infinite 

elastic  medium;  Phil.  Trans.  R.  S.  London  A,  841,  V.  242, 
63-100,  1949 

Newlands,  M.  The  disturbance  due  to  a  line  source  in  a  semi-infinite  elastic 

medium  with  a  single  surface  layer;  Phil.  Trans.  R.  S.  London 
245,  213-308,  1952 

Pekeris,  C.  L.  Theory  of  propagation  of  explosive  sounds  in  shallow  water; 

Mem.  Geol.  Soc.  Amert  27,  1948 

Sherman,  D.  I.  On  the  wave  propagation  in  a  liquid  layer  overlying  a  solid 

halfspace;  Pub.  Inst.  Seism.  Acad.  Sci.  USSR,  115,  1945 

Sommerfeld,  A.  Ueber  die  Ausbreitung  der  Wellen  in  der  drahtlosen 

Telegraphie;  Ann.  Phys.  28,  665-736,  1909 

"  "  Partial  Differential  Equations  in  Physics;  Acad.  Press,  New 

York,  1949 

Stone,  J.  L.  A  theoretical  analysis  of  acoustic  wave  modes  in  layered 

liquids;  Princeton  Univ.  Dept.  Electr.  Eng.  Tech.  Rep.  9, 

1953 

Whittaker,  E.  T.  and  Watson,  G.  N.  Modern  Analysis,  1947 


